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Abstract 


In this paper, the structural, electronic, and elastic properties of molybdenum disulfide (MoS2) 
in hexagonal and trigonal phases were investigated using first-principles calculations based on 
density functional theory (DFT) as implemented in the Quantum ESPRESSO package. Our 
findings from the structural and electronic properties of MoS2 show that bulk hexagonal and 
trigonal MoS phases have a calculated bandgap of 1.37 and 1.56 eV respectively, which is in 
good agreement with available theoretical and experimental results. The elastic properties show 
that both space groups of MoS, are brittle at zero pressure and they also agree with Born's 
mechanical stability condition. The overall results place our new optimized MoS:z as a potential 
candidate for optoelectronic application. 


Keywords: DFT; Transition metal dichalcogenides; anisotropy; chalcogen; elastic constants. 


I. INTRODUCTION 


1. metal dichalcogenides (TMDs) have piqued 
research interests globally due to their promising 
applications in the area of optoelectronics [1], spintronics [2], 
and nanoelectronics [3]. MoSz has been identified as one of 
the momentous 2D materials for the study of semiconductors 
and their application in physics and engineering [4], since its 
discovery. Graphene was the first thoroughly investigated 2D 
material [5], which offers distinct advantages because of its 
high electron mobility and thermal conductivity which made 
it suitable for electronic and mechanical strength applications 
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[6]. Graphene is a hexagonal carbon sheet with a two- 
dimensional structure and is immensely popular, but its zero 
bandgaps have sparked a search for 2D materials with 
semiconducting capabilities and applications [7-9]. 

TMDs are semiconductors materials of the general form 
MX2(M=Mo, W, etc.; X = S, Se, or Te) [10], the M is metal 
bonded to the non-metal called cuhalcogen X. Recently, 
revitalized interests of scientists and engineers in studying 
their structural, electrical, and elastic properties boosted to 
investigate its potential candidacy for the next generations of 
sustainability science and technology. TMDs exist in different 
phases with distinct space group numbers, with a lot of novels 
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that can also be used for many applications. Recent studies 
have been conducted on different phases of these materials 
and are found to be promising candidates for optoelectronic 
properties. For example, [11] studied the structural, electronic, 
and optical properties of WS2. Two bulk phases of WS2 with 
hexagonal phases and space group (P63/mmc) and 
rhombohedral phase (R3m) have been investigated and 
reported to be structurally stable both thermodynamically and 
dynamically. Reference [12] also investigated the electronic 
and optical properties of the hexagonal structure of the MoS2 
crystal for bulk and monolayer of MoS. A subsequent study 
by [13] used a self-consistent plane-wave method based on the 
DFT technique as implemented in the Quantum ESPRESSO 
using GGA-PW91, and their result demonstrates that the 
calculated lattice parameters overestimate the experimental 
values, which is a fundamental characteristic of standard GGA 
functional differing from experimental values [14, 15]. 
Reference [16] presented a study on the electronic, structural, 
and optical properties of graphene-like MoX2 using FP- 
LAPW+lo method in the WIEN2K simulation code. 
Similarly, [17] investigated the elastic properties of 
molybdenum disulfide using first-principles calculations, the 
authors calculated the elastic constant C;; of MoS: at a pressure 
ranging from 0 to 40 GPa, while [18] presented a study on the 
electronic and elastic properties of hexagonal layered crystals 
using van der Waals functional, with results indicating a 
monotonical increase in the elastic constant of MoSz and an 
increase in the anisotropies of MoS2 with an increase in 
pressure. 

Research on MoS2 has focused on the hexagonal space 

group's electronic, structural, and elastic capabilities, with 
little attention paid to the trigonal space group. However, 
earlier research has revealed that MoSz2 exists in different 
phases [19], with the graphene-like hexagonal bulk MoSz2 
being given more attention [20]. Furthermore, to the best of 
our knowledge, there is no comprehensive study of the 
structural, electronic, and elastic properties of a trigonal phase 
of bulk MoS2. As a result, research into the electrical and 
elastic characteristics is required. 
The present work is arranged as follows: An overview of the 
extensive computational framework used in this work is 
provided in Section 2. In Section 3, the results are illustrated 
together with physical explanations of the elastic, structural, 
and electrical characteristics of the materials under study. In 
Section 4, final remarks were made. 


II. COMPUTATIONAL FRAMEWORK 


In this work, the Materials Project (MP) database served as 
the primary source of the crystal structures [21]. A density 
functional theory (DFT) [22] based calculations have been 
applied using the generalized gradient approximation (GGA) 
[23] as implemented in the Quantum ESPRESSO package [24, 
25] employing the Perdew-Burke-Ernzerhof (PBE) variant for 
the exchange and correlation potentials. To improve the 
computation efficiency, the crystal structures of the two space 
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groups were fully relaxed at zero pressure until the forces on 
each atom were less than 0.01 eV/A, in which the optimized 
coordinates of MoS have been used for the calculations. A 
12 x 12 x 3 Monk Horst-Pack [26] of k-points was used for 
the sampling in the Brillouin zone. The spacing of the real 
space grid used to calculate the Hartree, exchange, and 
correlation contribution to the total energy and Hamiltonian 
was 400 Ry. A supplementary package to the quantum 
ESPRESSO thermo_pw [27] was also used in the calculation 
of the elastic properties using an energy cut-off of 60 Ry in 
conjunction with 12 x 12 x 3 for both hexagonal and the 
trigonal space group. 


III. RESULTS AND DISCUSSIONS 


A, Equilibrium Lattice Parameters 

Molybdenum disulfide has a hexagonal P63/ 
mmc [no194] with lattice parameters a = 3.19A andc = 
13.98 A and trigonal P3m1 [no164] with a lattice parameter 
of a = 3.19A and c = 24.88 A. The structural model of MoS2 
comprises Mo and S layers as shown in Fig. 1(a) and (b), for 
both hexagonal and trigonal space groups respectively [28]. 


Fig. 1 Structural model of bulk MoS; in (a) hexagonal and 
(b) trigonal phase 


The calculated structural properties of the bulk hexagonal 
and trigonal MoSz is presented in Table I and are compared 
with the available theoretical and experimental results. The 
bulk MoS: has two layers of Mo atoms; one layer is directly 
above the sulphur atoms of the other layer and vice-versa. 
After relaxation, the results of the optimized lattice parameters 
of MoS2 were compared with recently available calculated 
values. 

Table I. Calculated structural properties of bulk hexagonal and 
trigonal MoS2 


Bulk MoS; XC a(A) c(A) V(A‘) Bo(Gpa) By Ref. 

Hexagonal GGA 3.080 13.090 119.41 128.8 447 This study 
GGA 3.190 12.310 106.1 50.86 - [29] 
Exp 3.160 12.290 107.18 47.65 E [30] 
rvvl0 3.218% 12.448 [31] 

Trigonal GGA 3.45 21.470 194.51 76.6 4.56 This study 
GGA 357. 1838 - 2 Z [32] 
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To investigate the structural stability, the total energy of the 
two systems was fitted to volume using the Birch-Murnaghan 
equation of state, shown in (1) with the result of the fit in Table 
II and plotted in Fig. 2. 


Table II. Calculated bulk modulus, volume, pressure 
derivative, and total energy of MoS2 using the Murnaghan 
equation. 
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Fig. 2(a) the calculated energy against the volume of trigonal 
Mos: (b) pressure-volume graph of trigonal MoS2 


To determine the reliable theoretical values of the 
thermodynamical state of the trigonal phase, we use the total 
energy against the volume which is more reliable than the 
pressure against the calculated volume as can be seen in Fig. 
2(b), by applying the equation of state (EOS) at low cost. The 
calculated energy-volume fit fully demonstrates the trigonal 
space group of MoS> is thermodynamically stable as can be 
seen in Fig. 1(a). 


2 2 2 2 2 
E(V) =E,+ oe fe - i} Bit {i = i} {fe = +} (1) 


Where Ep is the ground state energy, Vp is the equilibrium 
volume and v is the deformed volume, By is the bulk modulus, 
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and Bg is the derivative of the bulk modulus. 


B. Electronic band structure and density of states 


To demonstrate the nature and the orbital contribution of 
the total density of state (TDOS) of the electronic energy band 
gaps of the two space groups MoS2, the knowledge of the total 
density of state and partial density of the state is critical. In 
this regard, the calculated bandgaps for the two phases of the 
material have been recorded alongside the experimental 
results in Table II. 


Table II. Calculated energy bandgap of hexagonal and 
trigonal bulk of MoS2 


Bulk MoS; xc bandgap (eV) Ref. 
Hexagonal GGA 1.37 This study 
1.23 [33] 
1.29 [34] 
rvvl0 0.85 [35] 
0.89 [36] 
Trigonal GGA 1.56 This study 


The electronic band structure and corresponding density of 
states or the hexagonal and trigonal bulk phases of MoS are 
shown in Fig. 3(a) and Fig. 4(a) respectively using GGA 
exchange-correlation, following the high symmetry 
points  — M — K —T, in the first Brillion zone, the partial 
and total density of state of hexagonal and trigonal were also 
shown in Fig. 3(b) and 4(b) respectively. The electronic band 
structure and density of states of the bulk hexagonal MoSz are 
divided into three energy states from the valence to the 
conduction bands spanning from —7.0 to +7.0 eV states 
respectively. For the hexagonal structure, the first group of 
bands in the electronic band structure and density of states 
around —7.00 to —1.00 eV is mainly due to the 3p and 4d 
states of S and Mo atoms [32]. The second group is Fermi 
energy in the energy range from —1.0 to 0.6 eV with no state 
contributing. In the third group above the Fermi energy, the 
energy range is 0.6 to 7.0 eV, the main contribution is due to 
Ad orbitals of Mo, and 3p of the S atom which is the main 
contribution, and it shows a strong hybridization. Similarly, 
the trigonal electronic band structure and density of states 
around —7.00 to —1.00 eV in the valence band are mainly due 
to the 3p and 4d state of S and Mo atoms contributing to the 
peaks in the TDOS [33]. The second group in the energy range 
from —1.00 to 0.75 eV within the Fermi energy state holds the 
calculated bandgaps of the trigonal bulk MoS) materials. In 
the third group above the Fermi energy, the energy range is 
between 0.75 eV and above in the conduction bands are 
mainly contributions from the 4d state of Mo and 3p energy 
state of S. The PDOS analysis demonstrates high contribution 
by 4d and 3p states of Mo and S respectively. 
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Fig. | Calculated (a) electronic energy bandgap and (b) density of states of hexagonal MoS» using GGA 
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Fig. 2 Calculated (a) electronic energy bandgap and (b) density of states of trigonal MoS2 using GGA-PBE 
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To further demonstrate the electronic band gap, the bands 
around the forbidden gap are relatively flat, as expected from 
the ‘d’ orbital. For bulk hexagonal MoS:, the valence band 
maximum is at I and the conduction band minimum is in 
between K and I direction of the hexagonal phase of MoS2 
revealing an indirect band-gap semiconductor, as can be seen 
in Fig. 3(a). For trigonal bulk MoSz, the valence band 
maximum is also at I’ and the conduction band minimum is at 
K high symmetry point revealing an indirect bandgap 
semiconductor as can also be seen in Fig. 4(a). Although our 
calculated band gaps for bulk MoSz are in good agreement 
with theoretical values obtained in [33], they underestimate 
the experimental band gap due to the inherent drawback of 
standard GGA functionals. Our calculated values of the band 
gap for hexagonal bulk MoS2 using GGA are 1.37 eV and 
1.52 eV for trigonal bulk MoS, which is in very good 
agreement with the results obtained in [12, 31-33]. In addition, 
it is found that the TDOS of hexagonal and trigonal phases are 
slightly similar as can be seen in Fig. 2(b) and 3(b), this is 
because both the two phases of MoS2 are mainly composed 
with Mo 4d, S 3p, and S 3s state to the TDOS. The added 
conduction band of hexagonal MoS: migrates from the high- 
energy region to the Fermi level, as a result, it promotes the 
electron transfer between the valence and the conduction 
band. This gives us a clear understanding of why the band gap 
of the hexagonal phase is smaller than the trigonal phase. 


C. Mechanical properties 


Several significant solid-state properties of crystals are 
connected to elastic properties, such as the equation of states, 
Debye temperature, mechanical properties, etc. Elastic 
modulus and elastic constants can provide a complete 
understanding of the structural behaviors of materials such as 
the stability of a material, hardness, ductility, brittleness, and 
other mechanical properties. In this Section, we focused on the 
crystal stability of the two space groups of MoS». To study and 
compute the elastic constants of these systems, the system 
responses to the applied stress need to be studied and analyzed 
from which all the elastic constants can be determined. In this 
work, we used Thermo-pw, auxiliary software in Quantum 
ESPRESSO code to calculate the responses of MoS» to 
applied stress. Molybdenum disulfide exists in hexagonal and 
trigonal space groups and both have five independent elastic 
constants Cy, Cy2, C13, C33, and C44, there exists in the 
calculations C6 [34] which is dependent on these five elastic 
constants, given by Cog = (Cy, — Cy2) / 2. Our calculated 
elastic constants for both hexagonal and trigonal are presented 
in Table 3 alongside available theoretical and experimental 
results [35]. The results obtained from the hexagonal space 
group are in total agreement with the experimental results in 
[36]. C72 and C73 are moderately close to the theoretical result 
obtained by [17]. To the best of our knowledge, there is no 
available result to compare with the trigonal space group 
MoSz. Under isotropic pressure, the mechanical stability of the 
two phases of both hexagonal and trigonal systems obeys the 
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Born mechanical stability criteria [37], as given by (1) to (3). 
The calculated elastic constants shown in (4) satisfy the Born 
stability criteria, suggesting that both phases of MoSz2 are 
mechanically stable at zero pressure. 


6 PU C. SOC, So (1) 
(Cy — Cy) > 0 (2) 
(Cy, + Cyz) > Cg3 — 2CZ3 > 0 (3) 


Table IV shows that both hexagonal and trigonal elastic 
constants are positive, indicating that the two compounds are 
stable. The ratios C33/C;; and C73/Ci2 show that atomic 
bonding along the x-direction is stronger than that along the z- 
direction for both hexagonal and trigonal atoms as can be seen 
in (3). Currently, to the best of our knowledge, there is no 
experimental or theoretical data on the trigonal space group 
for comparison. 


Table IV. Calculated elastic constants at zero pressure using 
GGA-PBE for bulk hexagonal and trigonal MoS2 


Phase XC Cr Cp Cp Cy Cy GyCn CyCn Ref 
Hexagonal 191 465 15 7 0.08  -0.11 This study 
249 3 37,30 —- : [17] 
238 +4 4 57 Oe : [35] 
EXP 238 -54 23 §2 19 - : [36] 
Trigonal 124 29 3 7 2 0.07 0.10 This study 


D. Elastic modulus and hardness 


The bulk modulus B together with the C4 can assess the 
machinability of the compounds via the machinability index 
H=B/C44, u provides the value that determines the degree to 
which a particular compound will be cut into different shapes. 
The higher value of yz the better the machinability. In this 
work, we have estimated the bulk modulus (8), the shear 
modulus (G), and the young modulus (Y) for the individual Cy 
for both hexagonal and trigonal structures respectively to 
determine their effect on the physical properties of the 
materials. Table V shows the calculated mechanical properties 
of the two structures. While the bulk modulus (B) describes 
the average atomic bond strength, the shear modulus (G) 
measures the resistance to restrict deformation upon shear 
stress. This demonstrates that there is a direct relationship 
between bulk modulus and shear modulus, implying that in the 
hunt for the hardness of a material, a large bulk modulus or 
shear modulus is required. The Pugh rule is used to estimate 
whether the MoS, will be brittle or ductile (B/G = 1.75) [39]. 
The proposed B/G ratio for determining material ductility and 
brittleness is around 1.75. If the B/G ratio is more than 1.75, 
the material is ductile; otherwise, it is brittle. In this work, the 
computed B/G values for hexagonal and trigonal space groups 
are 1.35 and 1.40, respectively, suggesting that both space 
groups are brittle at 0 GPa. To further demonstrate the nature 
of the bonding properties, we employ the service of the poison 
ratio, which tells information about the bonding of the system. 
For covalence, the poison ratio is usually less than 0.25 for all 
materials. Our findings in Table 4 show that the poison ratio 
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of the hexagonal and trigonal phases of MoS; is 0.20 and 0.21 
respectively. However, our results clearly show that both 
space groups of MoS; have covalent properties. Therefore, the 
covalent contribution for inter-atomic bonding is predicted for 
the two space groups of MoS. 


Table V. The calculated bulk moduli B in GPa, shear moduli 
G in GPa, Young modulus E£ in GPa, Pugh’s ratio G/B, and 
Poisson’s ratio at zero pressure. 


Phse B G Y BG v_ p=B/C4 Op —*Ref 
hexagonal 35 26 63 1.35 0.20 5.00 211.270 This study 
Hexagonal 51 47 107 —‘1,09 278.35 [17] 

Exp 53 

Trigonl 21 15 36 140 021 10,50 143,025 This study 


For the hexagonal or trigonal space group crystal structure of 

MoSz, the bulk, and the shear modulus are represented by (4) 

to (7). Equations (9) and (10) represent Young modulus (E) 

and Poisson’s ratio (v). 

Reuss’s approximation is given by (5) and (6) and Hill’s 

approximation is also given by (7) and (8). 

= [(C11+€13)C33-2C7] (5) 
2[3ByCaaCe6t{(C11+C12)C33+2C23}(CaatCaa)] 

_ (Cy +C13)C33-2C23 

C44 +Cq2 +2033 4+4C13 (6) 

G, = = (4611 + Cy2—-20334+4C,3412C44412C 56) (7) 

By = 5[2(Cia + Crz) + 4C,3—Ca3] (8) 


Young modulus E and Poisson's ratio can be deduced from the 


following relation: 
9BG 


ae 3B+G ) 
(10) 


y = 38226 
2(3B+G) 

The actual effective modulus of anisotropic polycrystalline 
crystals from the arithmetic means of these two values as 
approximated by Hill [40] the shear modulus and the bulk 
modulus are defined as: 

G == (Gp + G,) (11) 

B =~ (Br + By) (12) 
The Young modulus of a material is defined as the ratio of 
longitudinal stress and longitudinal strain, which is used to 
measure the stiffness of the solid material. A higher value of 
Young modulus indicates the stiffness of the material. From 
Table V, the hexagonal space group is stiffer than the trigonal 
space group. Another critical parameter that relates many 
physical properties of solids, such as their specific heat and 
melting temperatures is known as Debye temperature. The 
Debye temperature is obtained from the calculated elastic 


constants, which can be computed using (13). 
1 


Gr 


Br 


Sat = [eae |2 
Op — = *[=( M )| Vin (13) 
h is Planck's constant. 
-1 
1 2 3 
i= [aa] (4 
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G 
y= [exo (15) 
y= [eee (16) 


Where N,, is Avogadro's number, p is the density M, the 
molecular weight, V,, the average sound velocity, V, the 
longitudinal sound velocity andV, the transverse sound 
velocity. 


E. Anisotropy 

The orientation dependency of the elastic moduli or sound 
velocities is known as crystal elastic anisotropy [40]. The 
anisotropy can be defined by the three elastic independent 
shear constants. The hexagonal and trigonal shear anisotropy 


factors are given by (17) to (18) and the value of Re is K is also 
a 


computed. Table VI presents all the calculated anisotropy 
factors from independent elastic constants of hexagonal and 
trigonal phases [41]. 

1/6 (Cr+ C12 +2C33—4C13) 


A, = (17) 
C44 
2C44 
A, = ~—_ 18 
2 Cy -C12 (re) 
A; = ee (19) 
11 12 
Table VI. Anisotropy factors for bulk hexagonal and trigonal 
MoS. 
Phase Al A2 A3 K./Ka Ref 
Hexagonal 6.89 0.09 0.66 11.35 This study 
Trigonal 9.94 0.04 0.34 10.70 This study 


Table VI shows the hexagonal and trigonal phases of MoS2 
and the shear anisotropic factors calculated. A system with 
Anisotropy equal to | shows an isotropic property, while any 
value smaller or larger than 1 implies the material is 
anisotropic. The magnitude of the deviation from | is a 
measure of the degree of elastic anisotropy possessed by a 


crystal. Any value of Ke x. less than 1 reveals the 
a 


compressibility along the c-axis is smaller than that along the 
a-axis otherwise the compressibility along the c-axis is larger 


than that along the a-axis. Our calculated values for _ 7 K, for 
a 


both hexagonal and trigonal phases at zero pressure are 10.35 
and 10.07 respectively, the calculated directional bulk 
modulus of the two space groups suggest that their value are 
larger along the a-axis than on the c-axis, which indicate that 
the compressibility along the a-axis is smaller. 


IV. CONCLUSION 


The first-principle DFT calculation was used to calculate 
the structural, electronic, and elastic properties of bulk MoS2 
in hexagonal and trigonal phases as implemented in the 
Quantum ESPRESSO simulation code. Results obtained 
revealed that the hexagonal phase of MoS2 is much more 
stable than the trigonal space group. This study shows that 
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both space groups are semiconductors with an indirect 
bandgap of 1.37 and 1.56 eV for Hexagonal and trigonal 
phases respectively. The elastic properties demonstrated that 
both space groups satisfy Born's mechanical stability 
condition. The higher value of the Young modulus observed 
in the hexagonal phase indicates that the phase is stronger than 
the trigonal phase. Finally, the overall results place the new 
optimized MoS as a potential candidate for the next 
generation of sustainability science and technology to pave 
way for many technological applications. 
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